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ABSTRACT 

We describe global, 3D, time-dependent, non-radiative, general-relativistic, magnetohy- 
drodynamic simulations of accreting black holes (BHs). The simulations are designed to trans- 
port a large amount of magnetic flux to the center, more than the accreting gas can force into 
the BH. The excess magnetic flux remains outside the BH, impedes accretion, and leads to a 
magnetically arrested disc. We find powerful outflows. For a BH with spin parameter a = 0.5, 
the efficiency with which the accretion system generates outflowing energy in jets and winds 
is 7/ ^ 30%. For a = 0.99, we find rj ^ 140%, which means that more energy flows out of 
the BH than flows in. The only way this can happen is by extracting spin energy from the 
BH. Thus the a = 0.99 simulation represents an unambiguous demonstration, within an as- 
trophysically plausible scenario, of the extraction of net energy from a spinning BH via the 
Penrose-Blandford-Znajek mechanism. We suggest that magnetically arrested accretion might 
explain observations of active galactic nuclei with apparent t] ^ few x 100%. 

Key words: black hole physics — (magnetohydrodynamics) MHD — accretion, accretion 
discs — galaxies: jets — gamma-rays: bursts — methods: numerical 



1 INTRODUCTION 

Relativistic jets are 
holes (BHs). They 



a common 
are found 



feature of accreting black 

, ^ in both stellar-mass BHs 

jRemillard & McClintock l2006h and supermassive BHs in active 
galactic nuclei (AGN jTremaine et al.ll2002i) . Jets can be very pow- 
erful, with their energy output sometimes exceeding the Eddington 
limit of the BH. This suggests an efficient mechanism for their pro- 
duction. 

In seminal work, IPenrosd (1 19691) showed that a spinning BH 
has free energy that is, in principle, available to be tapped. This has 
led to the popular idea that the energy source b ehind relativistic jets 
is the rotational energy of the accreting BH. iBlandford & Znaie^ 
(Il97l hereafter IbZ7 7) came up with an astrophysical scenario in 
which this could be achieved. In their picture, magnetic field lines 
are kept confined around the BH by an accretion disc. The rotation 
of space-time near the BH twists these lines into helical magnetic 
springs which expand under their own pressure and accelerate any 
attached plasma. In the process, energy is extracted from the spin- 
ning BH and is transported out along the magnetic field, making a 
relativistic jet. The BZ mechanism is a promising idea since mag- 
netic fields are common in astrophysical accretion discs and so the 
requirements for this mechanism are easily met. 
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In the BZ mechanism, the rate at which rotational energy of 
the BH i s extracted - the BZ power Pbz - is given in Gaussia n-cgs 
units by (iBZVllTchekhovskov. Naravan & McKinnevll201Ql) 



Prz — 



47rc 



(1) 



where /c is a numerical constant whose value depends on the mag- 
netic field geometry (it is 0.053 for a split monopole geometry and 
0.044 for a parabolic geometry), Qh = ac/lru is the angular fre- 
quency of the BH horizon, Obh = (1/2) ^ {B'^ldAe^ is the mag- 
netic flux threading one hemisphere of the BH horizon (the integral 
is over all 6, (p at the BH horizon, and the factor of 1/2 converts 
it to one hemisphere), dAe^ = ^s^-gd6d^p is an area element in 
the 6-(f plane, and g is the determinant of the metric. Here a is 
the dimensionless BH spin parameter (sometimes also called a*), 
= ^g(l + Vl - a^) is the radius of the horizon, Vg - GM/c^ 
is the gravitational radius of the BH, and M is the BH mass. A 
simpler version of equation ([U with P oc was origina l ly de- 
rived by IbZVTI in the limit a «; 1. iTchekhovskov et all ilOld) 
showed that the modified form written here, with /(Qr) - 1. 
is accurate even for large spins up to a ^ 0.95, while for yet 
larger spins, they gave a more accurate 6th order approximation, 
/(Qh) ^ 1 + L38(QHr,/c)2 - 92(Qnrg/c)\ 

Using equation ([TJ, let us define the efficiency with which the 
BH generates jet power, t/bz, as the ratio of the time-average elec- 
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Figure 1. Shows results from the fiducial GRMHD simulation A0.99fc for a BH with spin parameter a = 0.99. The accreting gas in this simulation settles 
down to a magnetically arrested state of accretion. [Panels (a)-(d)]: The top and bottom rows show, respectively, equatorial (z = 0) and meridional (y = 0) 
snapshots of the flow, at the indicated times. Colour represents the logarithm of the fluid frame rest-mass density, log^opc^ (red shows high and blue low 
values; see colour bar), filled black circle shows BH horizon, and black lines show field lines in the image plane. [Panel (e)]: Time evolution of the rest-mass 
accretion rate, Mc^. The fluctuations are due to turbulent accretion and are normal. The long-term trends, which we show with a Gaussian smoothed (with 
width T = 1500r^/c) accretion rate, {M)jC^, are small (black dashed line). [Panel (f)]: Time evolution of the large-scale magnetic flux, 0bh, threading the BH 
horizon, normalized by {M)j . The magnetic flux continues to grow until t ^ 6000rg/c. Beyond this time, the flux saturates and the accretion is magnetically 
arrested. (Panels (c) and (d) are during this period). The large amplitude fluctuations are caused by quasi-periodic accumulation and escape of field line 
bundles in the vicinity of the BH. [Panel (g)] Time evolution of the energy outflow efficiency 77 (defined in eq.[5]and here normalized to {M)jC^). Note the large 
fluctuations in 77, which are well-correlated with corresponding fluctuations in 0bh- Dashed lines in panels (f) and (g) indicate time averaged values, <0|h^^^^ 
and (77), respectively. The average 77 is clearly greater than 100%, indicating that there is a net energy flow out of the BH. 



tromagnetic power that flows out of the BH, (Pbz>, to the time- 
average rate at which rest-mass energy flows into the BH, (M)c^, 



Vbz = 



(Pbz) 



X 100% 



Anc \ c 

\l/2 



^>/(Qh)x100%, (2) 



where 0bh = ^bh/ ((^>^Jc) is the dimensionless magnetic flux 
threading the BH and (...) is a time-average. Thus the efficiency 
with which a spinning BH can generate jet power depends on BH 
spin a via the angular frequency fin and on the dimensionless mag- 
netic flux 0BH- The strength of 0bh is very uncertain. 

It is generally agreed that 0bh is non-zero, since mag- 
netic flux is transported to the accreting BH by turbulent 
accretion. Howev er, the key elem e nts of this proce ss are 



not agreed upon jLubow et al.l Il994 ISpruit & Uzdenskvl l2005l : 



iRothstein & Lovel ace"2008': Beck with et al.ll2009l : ICadl201lh . This 
leads to a large uncertainty in the value of 77bz- Recent time- 
dependent general relativistic magnetohydrodynamic (GRMHD) 
numerical simulations have found a rather low efficiency, 77BZ ^ 



20%, even when the cen tral BH is nearly maximally spinning 
(McKinnev 2005; De Vill iers et all 12005: Hawley & Krolik 2006; 
liarkov & BaushevI 1201 ll) . With such a modest efficiency it is 
not clear that we are seeing energy extraction from the BH. 
The jet power could easily come from the accretion disc (see 
iGhosh & Abramowicz|[T997i : lLivio et aLlll999h . 

Observationally, there are indications that some AGN in 
the universe may have extremely efficient jets that require 
> 100% (lRawlings &Saundersl ll99ll: iGhisellini et aP I201Q| : 



iFernandes et all 1201 1[ iMcNamara et all 1201 ll : IPunslvl 1201 ih . A 
non- spinning BH usually has 77 < 10%, an d might under specia l 
circumstances have 77 ^ tens of percent (e.g. jNaravan et al.]|2003l) . 
However, a non-spinning BH can never give 77 > 100%, since this 
requires the system to produce more energy than the entire rest 
mass energy supplied by accretion. Values of 77 > 100% are possi- 
ble only by extracting energy from the spin of the BH. Thus, taken 
at face value, any robust observation of 77 > 100% in an AGN im- 
plies that the Penrose/BZ process must be operating. This raises the 
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following important question: Is it possible to show via a numerical 
simulation, an astrophysically plausible BH accretion scenario that 
gives jet efficiency 77 > 100%? To our knowledge, this has not been 
demonstrated with a GRMHD simulation. 

Here we describe numerical simulations in which we arrange 
our setup such that the accreting BH receives as much large- 
scale magnetic flux as can be pushed into the BH by accretion. 
In nonradiative MHD, the limiting flux is proportional to M. In 
fact, we supply more flux than this, so some of the flux remains 
outside the BH where it impedes t he accreting gas, le ading to a 
magnetically arrested disc" (MAD, Naravan et al.ll2003i, see also 



Bisnovatvi-Kogan & Ruzmaikinlll974 1 19761 : Igumenshchev et al.l 
2003h . The goal of the present simulations is to maximize 0bh and 
to make the jet efficiency as large as possible. As we show below, 
we do obtain larger efficiencies than reported in previous numerical 
experiments. Most interestingly, we find efficiencies greater than 
100% for a rapidly spinning BH. These experiments are the first 
demonstration of net energy extraction from spinning BHs via the 
Penrose/BZ mechanism in an astrophysically plausible setting. In 
^ we discuss our numerical method, the physics of MAD accre- 
tion and our problem setup, and in |3] we discuss the results and 
conclude. 



2 NUMERICAL METHOD AND PROBLEM SETUP 

We have carried out time-dependent simulations of BH accretion 
for two values of BH spin (see Table [B. We use the GRM HD code 
HARM dOammie et al.l20 03UMcKinn ey & Gammiel2004l) \yith re- 
cent i mprovements ([Mc Kinnev 2006: iTchekhovskov et al.ll2007l . 
I2OO9I : iMcKinnev & Blandford>(2009l) . The numerical method con- 
serves mass, angular momentum and energy to machine precision. 
We neglect radiative loss es, so the simulations corresp ond to an 
ADAF mode of accretion (iNaravan & McClintockll2008h . 

The simulations are carried out in spherical polar coordinates 
modified to concentrate resolution in the collimating polar jet and 
in the equatorial disc. We use Kerr-Schild horizon-penetrating co- 
ordinates. We place the inner radial boundary inside the BH outer 
horizon (but outside the inner horizon), which ensures that no sig- 
nals can propagate through the horizon to the BH exterior. The 
outer radial boundary is at r = lO^r^, which exceeds the light travel 
distance for the duration of the simulation. Thus both boundaries 
are causally disconnected. We use a logarithmically spaced radial 
grid, dr/(r-ro) - constant for r < rbr (see Table[T]for values of rbr), 
where we choose tq so that there are 9 grid cells between the inner 
radial boundary and the BH horizon. For r > rbr, the radial grid 
becomes progressively sparser, drir = 4(log r)^/^, with a smooth 
transition at rbr- At the poles, we use the standard reflecting bound- 
ary conditions, while in the azimuthal direction, we use periodic 
boundary conditions. In order to prevent the ^-extent of cells near 
the poles from limiting the time step, we smoothly deform the grid 
a few cells away from the pole so as to make it almost cylindrical 
near the BH horizon; this speeds up the simulations by a factor > 5. 

Numerical MHD schemes cannot handle vacuum. Therefore, 
whenever the fluid-frame rest-mass energy density, pc^, falls below 
a density floor PqootC^ = /?mag/^max, where p^ag is the magnetic 
pressure in the fluid frame, or when the internal energy density, Wg, 
falls below Wg^oor = O.lpfloor^^, we add mass or internal energy in 
the frame of a local zero angular momentum observer so as to make 
P = Pfloor or Ug = Wg^oor (McKinncy & Blandford 2009). The factor 
^max sets the maximum possible Lorentz factor of the jet outflow. To 
investigate the effect of this factor on jet efficiency, we have tried 



two values, ^max = 25, 250. There is little difference in the results. 
In any case, we track the amount of mass and internal energy added 
in each cell during the course of the simulation and we eliminate 
this contribution when calculating mass and energy fluxes. 

Model A0.99f (Table [B uses a resolution of 288 x 128 x 64 
along r-, 6-, and (f-, respectively, and a full azimuthal wedge, 
A(p = In. This setup results in a cell aspect ratio in the equato- 
rial region, 6r : rSO : r6(fi ^ 1 : 0.4 : 5. To check convergence 
with numerical resolution, at ^ = 14,674r^/c, well after the model 
reached steady- state, we dynamically increased the number of cells 
in the azimuthal direction by a factor of 2. We refer to this higher- 
resolution simulation as model A0.99fh and to A0.99f and A0.99fh 
combined as model A0.99fc. We also ran model AO. 99 with a 
smaller azimuthal wedge, A(p = n. We find that the time-averaged 
jet efficiencies of the four A0.99xx models agree to within statis- 
tical measurement uncertainty (Table [B, indicating that our results 
are converged with respect to azimuthal resolution and wedge size. 

Our fiducial model A0.99fc starts with a rapidly spinning BH 
(a - 0.99) at the center of an equilibrium hy drodynamic torus 
( Chakrabartilll985l : iDe VilHers & HawlevI l2003h . The inner edge 
of the torus is at = 15r^ and the pressure maximum is at 
^max = 34r^ (see Fig.[Tt). At r = rmax the initial torus has an aspect 
ratio h/r ^ 0.2 and fluid frame density p = 1 (in arbitrary units). 
The torus is seeded with a weak large-scale poloidal magnetic field 
(plasma/? = Pgas/Pmag > 100). This configuration is unstabl e to the 
magnetorotational instability (MRL lBalbus & HawlevI 1 99 ih which 
drives MHD turbulence and causes gas to accrete. The torus serves 
as a reservoir of mass and magnetic field for the accretion flow. 

Equation ([T]) shows that the BZ power is directly proportional 
to the square of the magnetic flux at the BH horizon, which is de- 
termined by the large-scale poloidal magnetic flux supplied to the 
BH by the accretion flow. The latter depends on the initial field 
configuration in the torus. Usually, the initial field is chosen to fol- 
low isodensity contours of the torus, e.g., the magnetic flux func- 
tion is taken as Oi(r, ^) = Cip^(r,6), where the constant factor Ci 
is tuned to achieve the desired minimum value of j3 in the torus, 
e.g., min/? = 100. The resulting poloidal magnetic field loop is 
centered at r = r^ax and contains a relatively small amount of mag- 
netic flux. If we wish to have an efficient jet, we need a torus with 
more magnetic flux, so that some of the flux remains outside the 
BH and leads to a MAD state of accretion jlgumenshchev et al] 
I2OO3I : iNaravan et al.ll2003l) . We achieve this in several steps. We 
consider a magnetic flux function, 0(r, 6) = r^p^(r, 6), and normal- 
ize the magnitude of the magnetic field at each point independently 
such that we have = constant everywhere in the torus. Using 
this field, we take the initial magnetic flux function as 02(r, ^) = 

^^le'-o IJ'-o" -^^^^^'^p' ^^^^ ^2 such that minfi = 100. This 
gives a poloidal field loop centered at r ^ 300r^. The loop has a 
much larger spatial size and more magnetic flux than the usual ini- 
tial field loop configuration considered in other studies. We main- 
tain a nearly uniform radial distribution of for r < SOOr^ which 
lets us resolve the fastest growing MRI wavelength with more than 
10 cells over a wide range of radii. 

Note that the above technique of starting with a large amount 
of poloidal magnetic flux in the torus is just a convenient trick to 
achieve a MAD state of accretion within the short time available in 
a numerical simulation. In our simulations magnetic flux is rapidly 
advected to the center from a distance r < lO^r^, whereas in na- 
ture we expect magnetic flux to be advected from a distant external 
medium at r > lO^r^ and to grow on the corresponding accretion 
time. The latter time is much too long to be currently simulated on 
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Table 1. Simulation details 



Name'^ 




n r%i 

'/ L /f^J 




Re^soliitinn (N v Nn ^ N ^ 


Smax 




^max/^g 


/? /r Pit r 




*run L'g/<-J 


*avg L'g/<-J 


A0.5 


0.5 


30 ±5 


n 


288 X 128 X 32 


25 


15 


34.475 


0.2 


200 


(0; 13,095) 


(10,300; 13,095) 


A0.99 


0.99 


145 ± 15 


n 


288 X 128 X 32 


25 


15 


34 


0.2 


200 


(0; 13,370) 


(6,000; 13,370) 


A0.99f 


0.99 


150 ± 10 


In 


288 X 128 X 64 


250 


15 


34 


0.2 


1000 


(0; 14,674) 


(7,000; 14,674) 


A0.99fh 


0.99 


135 ± 10 


In 


288 X 128 X 128 


250 


15 


34 


0.2 


1000 


(14,674; 30,500) 


(14,674; 30,500) 


A0.99fc 


0.99 


140+15 


In 


288 X 128 X {64, 128} 


250 


15 


34 


0.2 


1000 


(0; 30,500) 


(7,000; 30,500) 



^ Model A0.99fh is similar to A0.99f but with A^^ increased by a factor of two at t - I4,614rg/c. Model A0.99fc is comprised of models A0.99f and A0.99fh. 



a computer, hence the need to speed up the process in the simula- 
tions. Note that field advection is likely to be more efficient in a 

I 1 

thick , ADAF-like d isc rather than in a thin disc (e .g., Lubow et al. 
19941: ISpruit & Uz denskv 2005: Rothstein & Lov elace 2008; Cao 
20111) . hence we consider the present simulations to be a reason- 



able proxy for real ADAFs in nature. In any case, the key point 
is that, regardless of how a given system achieves a MAD state 
of accretion — whether it is by slow advection of field from large 
distances in a real system or through rapid advection of magnetic 
flux from short distances in our simulations — once the system has 
reached this state we expect its properties to be largely insensitive 
to its prior history. We thus believe the results obtained here are 
relevant to astrophysical objects with MADs. 



3 RESULTS 

Figure [T] shows results from the fiducial simulation A0.99fc for a 
rapidly spinning BH with a = 0.99. The top two panels (a) on the 
left show the initial torus, with purely poloidal magnetic field. The 
succeeding panels show how the accretion flow evolves. With in- 
creasing time, the MRI leads to MHD turbulence in the torus which 
causes the magnetized gas to accrete on the BH. In the process, 
magnetic flux is brought to the center and accumulates around the 
BH in an ordered bipolar configuration. These field lines are twisted 
into a helical shape as a result of space-time dragging by the spin- 
ning BH and they carry away energy along twin jets. 

The rate of accretion of rest mass, M(r), and rest mass energy, 
FM(r) = M(r)c^, at radius r are given by 



M(r) 



=-//' 



pu'dAe^ = FM(r)/c^ 



(3) 



where is the radial contravariant component of the 4-velocity, 
and the integral is over all 6, if at fixed r. The negative sign means 
that the flux is defined to be positive when rest mass flows into 
the BH. Figure [He) shows the rest mass energy flux into the BH, 
Fnirn) = M(rii)c^, as a function of time (the flux has been cor- 
rected for density floors, see Until a time t ~ 2000rg/c, the 
MRI is slowly building up inside the torus and there is no signifi- 
cant accretion. After this time, FM(ru) steadily grows until it satu- 
rates at ^ ~ 4000rg/c. Beyond this time, the accretion rate remains 
more or less steady at approximately 10 code units until the end of 
the simulation at ^ ~ SOOOOr^/c. The fluctuations seen in Fm are 
characteristic of turbulent accretion via the MRI. 

Figure \\lf) shows the time evolution of the dimensionless 
magnetic flux 0bh at the BH horizon. Since the accreting gas con- 
tinuously brings in new flux, 0bh continues to grow even after Fm 
saturates. However, there is a limit to how much flux the accre- 
tion disc can push into the BH. Hence, at ? ~ 6000r^/c, the flux 
on the BH saturates and after that remains roughly constant at a 
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Figure 2. Time-averaged fluxes of rest mass energy {Fm) and total energy 
{Fe) as a function of radius for models A0.5 (panel a) and A0.99fc (panel 
b). A positive flux means that the corresponding mass/energy flux is into the 
BH. For model A0.5 (a = 0.5), both Fm and Fe are positive, i.e., inward. 
The diff'erence (Fm - Fe) is the energy flux returned to infinity. The energy 
outflow efficiency (eq.[5} for this model is (t/) ^ 30%, i.e., 30% of the ac- 
creted rest mass energy is returned as outflow energy. For model A0.99fc, 
which has a rapidly rotating black hole (a = 0.99), (Fe) < 0, i.e., net energy 
flux is out of the BH. Correspondingly, {tj} ^ 140%. Greater than 100% effi- 
ciency is possible only by extracting net energy from the spinning BH. The 
gravitational mass of this BH decreases with time, though its irreducible 
mass increases. The effects of density and internal energy floors have been 
subtracted from the calculated fluxes. The subtraction is imperfect near the 
BH horizon, causing the slight non-constancy of fluxes at small r. 



value around 47. The correspondi ng dimensionless mag netization 
parameter T (iGammie et al.lll99^ : see |Penna"et d for defini- 

tion) is ^ 9.5 (much greater than 1), indicating that the flow near the 
BH is highly magnetized. Panel (b) shows that magnetic fields near 
the BH are so strong that they compress the inner accretion disc 
vertically and decrease its thickness. The accreting gas, of course, 
continues to bring even more flux, but this additional flux remains 
outside the BH. Panels (c) and (d) show what happens to the ex- 
cess flux. Even as the gas drags the magnetic field in, field bun- 
dles erupt outward (Ig umenshc hev 2008), leaving the time-average 
flux on the BH constant. For instance, two flux bundles are seen at 
X ~ +20rg in Figure [He) which originate in earlier eruption events. 
Other bundles are similarly seen in Figure [Hd). During each erup- 
tion, the mass accretion rate is partially suppressed, causing a dip 
in Mc^ (Fig.[T^); there is also a corresponding temporary dip in 0bh 
(Fig. [if). Note that, unlik e in 2D (axisymmetric) simulations (e.g., 
|Proga & Begelman"2003'). there is never a complete halt to the ac- 
cretion (Igumenshchev et al. 2003) and even during flux eruptions, 
accretion proceeds via spiral-like structures, as seen in FigurelHd). 

In analogy with Fm, let us define the rate of inward flow of 
total energy (as measured at infinity) as follows, 

FE(r)= f f TtdAe^, (4) 
Je J if 

where T^y is the stress-energy tensor. Figure [2] shows plots of 
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(FM(r)) and (^^(r)} vs r for the two simulations, A0.5 and A0.99fc. 
The fluxes have been averaged over time intervals (10300 - 
13095)r^/c and (7000 - 30500)r^/c, respectively, to reduce the ef- 
fect of fluctuations due to flux eruptions. The time intervals have 
been chosen to represent quasi-steady magnetically-arrested accre- 
tion. The calculated fluxes are very nearly constant out to r = 20M, 
indicating that both simulations have achieved steady state in their 
inner regions. Consider first the results for model A0.5 with a = 0.5 
(Fig.Eli). We find (Fm) ^ 12 units and (Fe) ^ 9 units. The dif- 
ference between these two fluxes represents the energy returned 
by the accretion flow to the external universe. Our simulations are 
non-radiative, with no energy lost via radiation. Hence the energy 
outflow is entirely in the form of jets and winds. 

We define energy outflow efficiency 77 as the energy return rate 
to infinity divided by the time-average rest mass accretion rate: 



(Fm) 



X 100%. 



(5) 



For model A0.5, the efficiency we obtain, (77) ^ 30%, is much 
larger than the maximum efficiencies seen in earlier simulations 
for this spin. The key difference is that, in our simulation, we max- 
imized the magnetic flux around the BH. This enables the system 
to produce a substantially more efficient outflow. 

In the more extreme model A0.99f with a = 0.99 (Fig.[2]3), we 
find (Fm) ~ 10 units and (Fe) ~ -4 units. The net energy flux in 
this simulation is out 6>/the BH, not into the BH, i.e., the outward 
energy flux via the Penrose/BZ mechanism overwhelms the entire 
mass energy flux flowing into the BH. Correspondingly, the effi- 
ciency is greater than 100%: (77) ^ 140%. Since the system steadily 
transports net total energy out to infinity, the gravitational mass of 
the BH decreases with time. Where does the energy come from? 
Not from the irreducible mass of the BH, which cannot decrease 
in classical GR. The energy comes from the free energy associ- 
ated with the spin of the BH. The BZ effect, which has efficiency 
T/Bz ^ 135% (eq. El with (02^> ^ 47^ from Fig. [if and k = 0.044), 
accounts for most of the extracted energy. 

Since greater than 100% efficiency has been a long- sought 
goal, we ran model A0.99fc for an unusually long time (t > 
30000r^/c). There is no indication that the large efficiency is a tem- 
porary fluctuation (see Fig.[T^). As a further check, we calculated 
efficiencies for each of the runs, A0.99, A0.99f, A0.99fh (Table[T]), 
to estimate the uncertainty in 77. We conclude that (77) ^ 140 +15% 
and that an outflow efficiency > 100% is achievable with a fairly 
realistic accretion scenario. We note, however, that by changing the 
initial setup, e.g., the geometry of the initial torus and the topol- 
ogy of the magnetic field, it might be possible to obtain even larger 
values of (77). This is an area for future investigation. 

Our outflows are in the form of twin collimated relativis- 
tic jets along the poles and less-collimated sub-relativistic winds 
jLovelace 1976; Blandford & Payne 1982). The former are mostly 
confined to streamlines that connect to the BH, while the latter 
emerge mostly from the inner regions of the accretion flow. The 
bulk of the outflow power is in the relativistic component. The en- 
ergy outflow efliciency shows considerable fluctuations with time 
(Fig.[T^), reaching values as large as 77 > 200% for prolonged peri- 
ods of time, with a long-term average value , (77) - 140+ 15%. This 
may explain sources with very efficie nt jets jMcNamara et al.l201ll : 
iFernandes et aDl201ll : |Punslvll201lh . The quasi-periodic nature of 
the fluctuations in 77 suggests magnetically-arrested accretion as a 
possible mechanis m to produce low-frequency QP Os in accreting 
stellar- mass BHs (iRemillard & McClintockl l2006l) and variability 
in AGN (iGhisellini et al.l201Ql) and GRB outflows JProga & Zhand 



l2006h . Additional studies are necessary to ensure the convergence 
of variability properties with numerical resolution. 

We conclude that rapidly spinning BHs embedded in 
magnetically-arrested accretion flows can produce efficient out- 
flows with (77) > 100%. Such flows could be relevant for under- 
standing astrophysical systems with extremely eflicient jets. The 
fiducial model A0.99fc presented here, which is designed to mimic 
magnetically arrested systems in nature, has a net energy flux away 
from the BH and demonstrates that net extraction of energy out of 
an accreting BH is viable via the Penrose/BZ effect. 
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